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Abstract 

We study theoretically the photon-assisted tunneling through a quantum dot 
weakly coupled to Luttinger liquids (LL) leads, and find that the zero bias 
dc conductance is strongly affected by the interactions in the LL leads. In 
comparison with the system with Fermi liquid (FL) leads, the sideband peaks 
of the dc conductance become blurring for ^ < g < 1, and finally merge 
into the central peak for g < ^ {g is the interaction parameter in the LL 
leads) . The sidebands are suppressed for LL leads with Coulomb interactions 
strong enough, and the conductance always appears as a single peak for any 
strength and frequency of the external time-dependent field. Furthermore, 
the quenching effect of the central peak for the FL case does not exist for 
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The photon-assisted tunnehng (PAT) through a quantum dot weakly coupled to two 
Fermi liquid leads (FL-QD-FL) has been studied intensively in the last decade. 1^-0 When 
two time-dependent fields with the same form V cos{Qt) are applied on the leads symmet- 
rically, a series of additional peaks (called sidebands) occurs in the curve of the zero bias 
dc conductance vs. the gate voltage. These peaks locate at energies 6E = nhQ away from 
the original resonant level. The nth sideband will be quenched as a equal to zeros of the 
Bessel function J„, @] because the height of the nth peak is proportional to (a) (where 
a = 1^ is a dimensionless variable). Physically, all these features of PAT can be ascribed 
to characters of the electron occupations in the leads, or equivalently the characters of the 
density of states and the Fermi distribution of the electrons in the leads. 

One dimensional interacting electron systems behave as Luttinger liquid (LL), for 
that theory predicts the power-law behavior of the density of states (DOS) of the electrons. 
01 What will happen for the sidebands features if the two FL leads are replaced by two LL 
leads? We expect that these features may be strongly affected by the specific characters of 
the DOS in LL leads. With this idea in mind, we study in this letter the PAT through a 
QD weakly coupled to two LL leads, i.e., a LL-QD-LL system. 

In contrary with the FL-QD-FL system, we find that the sideband peaks of the dc 
conductance become blurring for ^ < g < 1, and finally merge into the central peak for 
g <\- The sidebands are suppressed for LL leads with Coulomb interactions strong enough, 
and the conductance appears as a single peak for any strength and frequency of the external 
time-dependent field. We call this the submergence of the sidebands. The contributions of 
each PAT process to the corresponding sidebands are diminished by the suppression of the 
DOS near Fermi energy, and shared due to the unusual DOS of the LL leads. Furthermore, 
in strong interaction case (f? < |) , the central peak consists not only of the direct tunneling 
process but also many other type PAT processes, so that it does not exist a single parameter 
a by which the central peak can be quenched. 

We use the nonequilibrium-Green-function technique to calculate the average current 
through the quantum dot. ^ Under the adiabatic approximation, the electron occupations 
in each lead remains unchanged. Thus the exact time-dependent Green functions in the 
uncoupled leads are 

9^t^itu t2) = exp(-z r Ai/^(t)rft)^7;;^';/(ti, t2), (1) 

Jt2 

where fi'o L/i?'^(^i' ^2) are the Green functions of the uncoupled leads without the time- 
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dependent field, and Ai//j(t) = eVL/RCOs{Qt). We take the units of h = 1 hereafter. 

The Hamiltonian of the system can be split into three parts: H = Hieads + Hd + Ht, 
where Hieads = Hl + Hr denotes the left and right LL leads, Hd = ed^d is the Hamiltonian 
of the QD, with {d^,d} the creation/annihilation operators of the electron in the dot and 
only one energy level is considered. Ht is the tunneling Hamiltonian and can be written as: 

Ht= {txd+^x + h.c), (2) 

\=L,R 

in which {ipx^ V'a}(A = L/R) are the Fermi operators at the end points of the left/right lead; 
{t\} is the tunneling constant. The spin indices have been suppressed since they are not 
important in the following discussion. 

Using nonequilibrium-Green-function technique, the current through the quantum dot 
can be formulated as: 

J(t) = JL{t) = -2eRe J dU[tl{G'^ {t,U) g< {h,t) + G< {t,h) g'l{h,t))tL] , (3) 

where we have defined G^ (t,ti) = ~i9(t — ti) {{d(t) ,d'^ and all the other Green 

functions with the standard forms. 

Supposing J(t) = -'n exp(mf2t), then 

J„ = -2e|ti|^Re^ / ^ (G^^^.o H ^L,n,„^ H + G< ,o M J^) , (4) 

in which the functions F„^m(u;) = Fn~m{l^ + rnVL), and Fn{uj) is the double Fourier transfor- 
mation of function F{t,t'): [Tl[] 

F{t,t) = ^exp(mfit') / pe-'^(''')F^{u). (5) 

27r 

The retarded Green function of the QD in equation (4), G^ it,ti), can be calculated by 
using Dyson equation 

G^ (t, tl) = (t, h) + J J drdr'g^ (t, r) (r, r') G'' [r m) , (6) 

where {t,t'^ = J2\=L,R\'tx\'^ 9\'^ {j'^'^'^ is the irreducible self-energy, g^ {t,ti) is the 
free retarded Green function of isolated dot. Notice that the Green functions of the QD 
depend on two time variables, not the time difference; so one should take the double Fourier 
transform, and obtain: 
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h,h 

Neglecting the off-diagonal elements of /jl"-^)) we have 

Qr _ ^"."t 

in which J; (a) is the Bessel function, e is the energy level in the dot and Ql/r 
lesser Green function can be calculated using Keldysh formalism: 



(7) 



(8) 



The 



(9) 



Substituting equation (8) and (9) into equation (4), we can obtain the average current 



as: 



{J) = e\tLtR\: 



duj (Em Ji{aL)9oA^ - "^fi))(Em Jl,{(^ r) 9^r{^ - mV)) -{L^ R) 
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(10) 



Equation (10) is the main result of our work. For the model discussed in this letter, the 



Green functions are taken from ||T2 

goj (^) 



T 



(11) 



90R l^J 



±z- 



T 



eH) 



where T is the temperature, and H the dc bias voltage. ■jL/Ri^) is defined as 'jL/Ri^) 



L/R 



2n\t 



L/R I 



(ir) ^^^^^ F(i/2gL/fl+»£/27rr)| ^ j^^^^ Gamma function, g^/ji is the inter- 

^ ^ ^[^/gL/R ) 

action parameters characterizing the left /right LL liquids, F^/^ describes the effective level 
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broadening of the dot, proportional to 



L/R 



and A is the high-energy cutoff or a band 



width. The Hamiltonian can describe a quantum dot coupled to the end points of either 
two half-infinite quantum wires or two edge states of fractional quantum Hall liquids with 
filling factor u = 2m+i ~ 9- '^^^ average current formula (Eq. (10)) covers a wide range of 
interaction in the leads, from noninteracting limit {g = l)to strong interacting LL (small g) . 

What we are interested in is the dc conductance at zero bias. We consider the symmetry 
case, namely = tji = t, gi = gR = g, and ai = ctR = a. If the external fields are not very 

\ 1/9-1 



strong, equation (10) can be further simplified. Under the condition ^ 



In ^ « 
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1(T < 



or 



1/9-2 



< i(r > 



[p!0| we obtain the zero bias conductance after a 
contour integral (all the subscript L/R are omitted): 

,2 (Em ^m(a)^(r(£ - "^^))£(Em Jl{a)9o{e - mQ)) - (<^>) 



G 



t 



(12) 



In the following, we use current formula (12) to study the zero bias conductance vs. 
the energy level of the QD, G/Gmax vs. ^, numerically. We consider the whole range of 
interaction of the leads, from noninteracting limit {g = 1) to strong interaction case {g < i), 
with much attention on the latter. Fig.l presents the results for p = 5, a = ^ = 1 and 
g = 1.0,0.8,0.5,0.2. The sideband peaks appears clearly for g = 1; become blurring for 
g = 0.8; merge into a single central peak for g = 0.5; and finally become narrower for 
9 = 0.2. 

In order to study the quenching behavior, we chose a equal to the smallest positive root 
of the zero order Bessel function, Jq. For the noninteracting limit case {g = 1), the central 
peak is completely quenched as expected. However, with the increase of the interaction or 
the decrease of g, the situation changes significantly. For g = 0.8, the relative heights of the 
sideband peaks are almost unchanged, but the valleys lifted. Both are shown in Fig. 2 (a). 
Most surprisingly, the quenching effect disappears completely for the strong interaction case 
(shown in Fig.2(b)): a broader peak with a flat plateau for g = 0.5; and a narrower single 
central peak for g = 0.2. 

To further understand the PAT in LL-QD-LL for strong interaction case, we calculate 
C/Gmax vs. j; for g = 0.2, with different strengths and frequencies of the external fields, 
shown in Fig. 3. One can clearly see that all sideband peaks disappear. And the width of the 
single central peak becomes broader with the increase of the strength or the frequency of the 
fields, originated from the contributions of the multi-photon-assisted tunneling processes. 

All features mentioned above can be ascribed to the power-law behavior of the DOS and 
the suppression of the electron tunneling near the Fermi energies of LL leads. In order to 
understand the underlying physics, we label (m, n) to indicate two pairs of processes: (l)the 
electron tunneling through the dot with m photons absorbed in the left and n photons 
emitted in the right lead, and its reversal process; (2) same processes as in (1) except m 
^ n. As an analogy of the (m, n) , one can consider the tunneling currents under two 
different finite bias voltages (mf2 + eVf,) — nQ and {nQ + eVf,) — mQ, without the time- 
dependent fields. The difference of these two net currents is approximately equal to the net 
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current of (m, n), except the effect due to difference of tfie width of the energy level in QD. 
In the following, we shall discuss the PAT quahtatively with the help of this analogy. 

It is well-known that for FL-QD-FL, the current under finite bias change abruptly at 
the Fermi energies in the two leads due to the abrupt change of the electron occupations 
near Fermi energies in the leads; thus two sharp peaks at £ = mfl and nfl obtained from 
the difference of the currents under two finite bias voltages, provided Vb is small. Similarly, 
the conductance from (m, n) has the same peaks at e = mQ and nQ. Conductance at other 
energies is strongly suppressed due to the uniform distribution of the electrons. 

However, for LL-QD-LL, the DOS of the leads behaves as p (a;) oc |a;|^^^~^ as T ^ 0, the 
so-called power-law behavior of the DOS, so p' {lo) oc . For the weak interaction ( 

g > \ p' {uj) — > oo as a; — > 0, thus there is also an abrupt change in the electron occupation 
near Fermi energy; two peaks at £ = nQ and mQ still remain, but the valley between them is 
hfted due to the nonuniform of the DOS. The contributions from all different (m, n) processes 
give the total conductance with sideband peaks at £ = ±Q, ±2Q..., and the valleys are hfted 
with the decrease of consistent with the curves for g = 1, 0.8 in Fig.l and Fig.2(a) . 

For the strong interaction case ( (? < | ), p' (uj) ^ as ^ 0, thus no abrupt change 
in the electron occupation exists, and the difference of the currents under two finite bias 
voltages rises slowly at the outset and finally appears as a single peak centered at ^^^Q 
approximately, so that the (m, n) processes do not contribute significantly at £ = mil and 
nJl, instead, a peak at £ = ^^^Q is produced. Furthermore, our calculation shows that 
processes such as (m, m), always producing a single peak at £ = mQ for any values of g, 
are not important, due to the fact that the direct tunneling near Fermi energy is strongly 
suppressed for strong interaction case. In Fig. 4, we present the results from processes (1, 1) 
and (2,2), showing a contribution almost 6 orders of magnitude smaller than the main 
processes. Moreover, one can see from Fig. 4 that the main contributions to the conductance 
are originated from processes (m, —n), with m,n = 0, 1, 2..., especially from that of m = n, 
which make significant contribution at £ = but little at the sideband energies. That is why 
for the strong interaction case only one single peak at £ = exists, while all the sidebands 
are suppressed completely, as shown in Fig.3. 

Now let us explain why the quenching of the central peak disappears for strong interaction 
case. Notice that the main contributions of the central peak come mainly from processes 
(0, m) with m e Z for weak interaction (^f > 1/2); but (m, — n) with m ~ n for strong 
interaction case ( g < 1/2). Given a the zeros of Jq, all the processes like (0, m) become 
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zero and contribute nothing to the central peak for weak interaction case, but the processes 
in the strong interaction case are not affected much. Thus the central peak is quenched 
as usual in the former, but will be lack in the latter (see Fig.2). Furthermore, since the 
conductance at e = come from many PAT processes, it does not exist a single parameter a 
by which all the processes,having main contributions to the central peak,can be neglected; 
thus there is no quenching of the central peak in case of LL leads with g < \- 

Formula (10) can also be applied to the case with asymmetric external fields. If only one 
lead is applied with an ac field, we find that the photon-electron pumping effect still occurs 
for LL case. However, the negative current becomes smaller and smaller with the decrease 
of g and the maximum of the current shifts towards the negative e (not shown here). 

In conclusion, we have studied the PAT in the system of LL-QD-LL for the whole range 
of interaction of the LL leads and the crossover behavior. We find that the zero bias conduc- 
tance is strongly renormalized by the interaction. In particular, for the strong interaction 
case( (7 < I ), all sidebands merge into a single peak (we called it submergence of the side- 
bands), and the quenching effect also disappears completely. These novel features of PAT 
have been ascribed to the power-law behavior of the DOS and the suppression of the elec- 
tron tunneling near the Fermi energies of LL leads. Since recent experiments have confirmed 
that the carbon nanotubes (CNT) do have some properties of the LL, [0,0] especially the 
power-law feature of the DOS, we suggest a setup consisted of CNT-QD-CNT to study the 
PAT of the system. We hope our theoretical predictions can be tested experimentally. 

This work was supported by NSFC under grant No. 10074001. One of the authors 
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scopic Physics in Peking University. 
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FIGURE CAPTIONS 

Fig. 1 G/Grnax VS. e/T for g = 1.0,0.8,0.5,0.2. G is the zero bias conductance, Gmax is the 
maximum of G, e is the energy level in the QD, T is the temperature, and g is the 
interaction parameters of the leads. The parameters of the time-dependent fields are: 
n/T = 5 and a = ^ = 1. 

Fig. 2 G/Grn!,x vs. e/T for (a) g = 1,0.8; (b) g = 0.5,0.2, with n/T = 10 and a the first 
positive root of Jo- 
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Fig. 3 G/Gmax vs. e/T for g — 0.2. The parameters are chosen in (a) a — 0,1,3,5 and 
Q/T = 5; (b) q; = 1 and Q/T = 0, 5, 10, 15. 

Fig. 4 The contributions of some {m,n) processes for g = 0.2, fl/T = 10, and a the first 
positive root of Jq. Gmax is the maximum value of the total conductance. 
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